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Abstract 

The determination of the elements of the S-matrix within the framework of 
time-dependent density-functional theory (TDDFT) has remained a widely open question. 



^ . We explore two different methods to calculate state-to-state transition probabilities. The 

| 7 | I first method closely follows the extraction of the S-matrix from the time-dependent 

^ ' Hartree-Fock approximation. This method suffers from cross-channel correlations resulting 

in oscillating transition probabilities in the asymptotic channels. An alternative method is 
I proposed which corresponds to an implicit functional in the time-dependent density. It 

^ gives rise to stable and accurate transition probabilities. An exactly solvable two-electron 

system serves as benchmark for a quantitative test. 

As a matter of principle, time-dependent density functional theory ^ provides a highly 
efficient method to solve the time- dependent quantum many-body problem. It yields directly 
■ the time-dependent one-particle density n(r, t) of the many-body system. All physical ob- 
Q^, servables of the quantum system can, in principle, be determined from the density. In 



practice there a two essential ingredients to a TDDFT calculation. First an approximation 
to the time- dependent exchange-correlation potential \4.c[n](r, t) has to be found which via 
the non-interacting Kohn-Sham system determines the evolution of the density. The second 
ingredient are functionals that allow the extraction of physical observables from the density. 
For some of the observables such as the ground-state energy extraction is straight forward 
within ground-state density functional theory . Excited-state spectra have been obtained 
from linear-response functionals 3, |^. Beyond linear response, the time-dependent dipole 
moment which governs the emission of high-harmonic radiation can be directly determined 
from n{f,t). Ionization probabilities can be approximately extracted by identifying the in- 
tegrated density beyond a certain critical distance from the bound system with the flux of 
ionized particles ^,0]. However, in general, on the most fundamental level, state-to-state 
transition probabilities contain the full information on the response of a many-body system 
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to an external perturbation. One example are bound-bound transition amplitudes rec[uired, 
e. g. in coherent control calculations of laser-matter interactions within TDDFT [7], cur- 
rently a hot topic since atto-second laser pulses allow the control of the electron dynamics. 
This poses one fundamental question: How can state-to-state transition probabilities be 
extracted from TDDFT? 

The ultimate goal of the study of the in general non-linear response of the many-body sys- 
tem to a time-dependent perturbation is the determination of the state-to-state transition 
amplitude 

S,f=\im{xf\Uit,-t)\x^), (0.1) 

I— »oo 

where \Xi,f) are the initial (final) channel states of the system prior to (z) and after (/) the 
perturbation, and U{ti,t2) is the time evolution operator of the system. The challenge is, 
thus, to construct a functional Sij[n] that allows to extract Sif from TDDFT. The present 
paper addresses methods to extract transition probabilities between discrete states of the 
many-body system. As point of reference, we investigate first the evaluation of eq. (jO.lll 
employing Kohn-Sham orbitals in close analogy to the time-dependent Hartree-Fock (TDHF) 

nnnn 

method (see [a, ISli llOl lUl and references therein). This method would be the equivalent of 
the at least zeroth order of a time-dependent many-body perturbation theory (S'-matrix 
theory) in which the time-dependent non-interacting Kohn-Sham system is considered as 
the unperturbed system. This method involves three steps of approximations: The initial 
state, the final state and the many-body propagator are approximated by their TDDFT 
equivalents. We encounter similar conceptual problems ("cross-channel correlations") as 
TDHF does. We then formulate a novel functional that allows the determination of Sij[n] 
which is shown to be free of these deficiencies. We test the method with the help of an 
exactly solvable two-electron model. 

We consider an interacting A^-electron system of Hamiltonian Hq with stationary eigenstates 
Xi.f which is subject to a perturbation V{t) which is switched on at time t = and switched 
off at time t = t. The initial state of the system \xi) is assumed to be the ground state and 
evolves according to the time-dependent many-body Schrodinger equation (in a.u.) 

z^^im) = [Ho + vmm) , i^(o)) = \x^) ■ m) 

The state-to-state transition amplitude (or S-matrix) from the initial state \xi) to a final 
state \xf) is defined by the overlap of the propagated state |^l/(t)) with eigenstates \xf) of 
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the unperturbed system 

S,j=\im{xf\^it)). (0.3) 

t — ^oo 

For later reference we note that the time evolution of the projection amplitude for t > r is 
given in terms of the eigenenergies of the asymptotic final states, £/, by 

{Xf\^it)) = expi-ze,it-T)){xf\^ir)) . (0.4) 

Since the perturbation vanishes for t > t, the state-to-state transition probability is given 

byP./ = |5./P = Kx/l^(r))p. 

Within TDDFT, the time-dependent density is represented through the time-dependent 
Kohn-Sham spin-orbitals ^fjj{f,t) as 

n(f,t) = J2 ^-(^'^) =J2Y1 ' (0-5) 

where Ncr denotes the number of electrons of spin cr. The one-particle spin-orbitals $o-j (^, t) 
evolve according to the time- dependent Kohn-Sham equation governed by the one-particle 
Kohn-Sham Hamiltonian 

+Vh [n] (f , t) + Kc [n^ , n^] (r , t) (0.6) 

which includes the external one-particle potentials, the Hartree potential and the exchange- 
correlation potential. The initial states |$o-,j(0)) = |$<j,j) are the occupied Kohn-Sham 
orbitals of stationary ground state density functional theory (DFT). Although Kohn-Sham 
orbitals have, a priori, no physical meaning as single-particle quantum states, the Slater 
determinant of Kohn-Sham orbitals, |\['™^^'^) := i, .., <l>| 7V|5 •■•5 '^'i.Af;), where 

A denotes the operator for anti-symmetrization, may be interpreted as zeroth-order ap- 
proximation to the many-body wavefunction in terms of coupling-constant perturbation 
theory 12. |l3|- It is therefore tempting to determine, in analogy to the TDHF approx- 
imation [8, U an approximate S-matrix as the projection amplitude Sij{t) ~ 
(X/|^™^^'^(t)). A dehcate quest ion arises at this point: Which are the appropriate chan- 
nel states Xij to project on? For the initial state, stationarity of the propagation of the 
system in the limit of a vanishing external perturbation {V{f, t) = 0) mandates that Xi is 
a Kohn-Sham Slater determinant of the occupied ground-state orbital. No such restriction 
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is imposed on when the evolution is calculated by forward-propagation. The simplest 
choice for channel states \xf) are Kohn-Sham Slater determinants built up from occupied 
and virtual Kohn-Sham orbitals |$ctj) of the ground state DFT problem, i.e. 

Sif^ lim(xr^|vl>™^^(t)). (0.7) 

Reliable transition probabilities can only be expected if both, time-dependent and station- 
ary Kohn-Sham Slater determinants are good approximations for the time-dependent and 
stationary many-body wavef unctions, respectively. Excited states often show a higher de- 
gree of correlation and a single Kohn-Sham Slater determinant is no longer a satisfactory 
approximation. In those cases, more elaborate final-state channel functions are needed. Al- 
ternatively, confi gura tion- interaction (CI) or multi-configuration Hartree-Fock (MCHF) can 
be employed jl^. llS^ . Within TDDFT, linear response theory allows the calculation of the 
excitation spectrum by including particle-hole excitations 0, As a by-product, improved 
excited-states, i.e. the particle-hole reduced density matrix, are generated in terms of an 
expansion in single-particle excitations of Kohn-Sham orbitals. In the present case of a 
two-electron system this approach should give an improved wavefunction compared to the 
initial single Kohn-Sham Slater determinant 16]. As discussed below, also this approach 
suffers from the same short-comings as does eq. ()0.7p . 

We therefore introduce a new functional which depends only on the time-dependent den- 
sity. For simplicity and in line with the present model system, the derivation is presented 
for two-electron systems. Generalizations to arbitrary many-electron systems are straight 
forward. Starting point is the expansion of the exact time-dependent wavefunction in terms 
of a complete set of final-state wavefunctions 

^(fi, f2, t) = ^(x/|^(t))x/(r-'i, fa) . (0.8) 
/ 

Using the stationary one-particle reduced density matrix 

P/^j(r) = 2 J df^2 x}'ir,f2)xf{r,f2) (0.9) 
and the time-dependent transition density matrix defined by 

Tf,jit) = {xM^it)y{xf\^it)) (0.10) 

the exact time-dependent density is given by 

n{f,t) = J2TfAt)P%irl = Tr [r(t)p«(r)] . (0.11) 
fj' 
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FIG. 1: Comparison of occupation probabilities of the ground state (figure a), first excited state 
I'^rei — 0) (figure b) and second excited state \N(~m — '^^^rei — 0) of exact calculation 
(full line) and TDDFT calculation (dashed line) for a confining frequency uj = 0.25. TDDFT 
occupation probabilities are obtained using the approximate S-matrix eq. HU.7() , final channel states 
are Slater determinants of Kohn-Sham orbitals |0, 0) ,[0, 1) and |0, 2). Parameters of the laser pulse: 
Fq = 0.07, UJL = 0.1839, r = 168 

The transition density matrix and , in particular, transition probabilities I'S'ijp can be 
directly determined by inversion of eq. ()0.11|) . Our primary interest lies in the diagonal 
elements (/' = /) of the transition density matrix Tf'j{t) —>■ S* j:,{t)Sij{t) at times after the 
switch-off of the external perturbation t > t. In this case the inversion problem of dimension 
Np X Np {Np: dimension of truncated final-state space considered) can be drastically sim- 
plified. Using eq. ()0.4|) . for non-degenerate final states (e/ 7^ Sf), the transition probabilities 
Tff can be extracted from a time-average over an interval {t — T)\ef> — Sf \ ^ 27r, 



n(r, t) := 




(0.12) 



leading to the read-out functional 

limn{f,t) = J2pf}mSM' . (0.13) 

^00 ^ 
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Unlike eq. ()U.lip . eq. ()0.13|) requires only an A^ir-dimensional inversion. In practice, the 
application of eq. ()0.13p requires evaluating the final state densities pf\{^ at Np distinct 
points fj j = l,...,Nf, so that the matrix Rfj := p^pj-{fj) does not become near-singular 
and remains invertible. The state-to-state transition probabilities then become 

|5,;r= limX]n(r^,t)i?7j /=l,..,iV/. (0.14) 

We have tested the functional (eqs. ()0.7|) and (10.141) f or an exactly solvable system of 
two electrons confined to a harmonic quantum dot In its present ID version, 

the electron-electron interaction must be replaced by a soft Coulomb potential The 
Hamiltonian of the system is given by 



where Xi and Pi are the coordinates and momenta of electron i {i = 1,2). The softening 
parameter b is set to b = 0.55. Similar systems to describe helium in one dimension were 



studied in the past 



m 



19|. The laser field F{t) with driving frequency ul = 0.1839 



and a peak field amplitude Fq = 0.07 is treated in dipole approximation. A pulse length 
r = 168 was chosen with a two cycle turn on, a two cycle fiat top and a two cycle turn 
off. Introducing center of mass (c.o.m.) coordinates R = xi + X2 and relative coordinates 
r = xi — X2 the Schrodinger equation ()0.2|) can be separated, since the Hamiltonian of eq. 
()0.15|) splits into H = H^ei + Hcm(t)- The eigenstates of the unperturbed system are char- 
acterized by the set of quantum numbers {Ncm, ^rei)) the number of nodes of the c.o.m. and 
relative wavefunction. Initial state is the spin-singlet ground state \Ncm = 0,nrei = 0). The 
time dependence of the total Hamiltonian is confined to the c.o.m. Hamiltonian. The exact 
time-dependent wavefunction therefore separates into a time-dependent c.o.m. and a time- 
independent relative part, \I'(r, i?,t) = g{r)h{R,t). Since the system starts out from the 
ground state, h{R, t) represents a coherent state of the one-dimensional harmonic oscillator 



namics of the density of this model system is 
20|,|2l|, i.e. the density is rigidly shifted without 



driven by an external electric field. The c 
subject to the harmonic potential theorem 
any distortion. 

In the case of a two-electron spin-singlet system evolving from the ground state, the time- 
dependent Kohn-Sham scheme consists in solving the Kohn-Sham equation for one doubly- 
occupied Kohn-Sham orbital $(x,t). For spin-unpolarized two-particle systems the exact 
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FIG. 2: Comparison of exact occupation probabilities (solid line) and transition probabilities ob- 
tained by inversion of eq. (|U.14|) (dashed line) for the ground state (figure a), first excited state 
(figure b) and second excited state (figure c). 

Vxc{t) can be constructed using tlie exact density derived from the Schrodinger equation and 
inverting the Kohn-Sham equation 2^. Since the harmonic two-electron quantum dot 
satisfies the harmonic potential theorem the time- dependence of the exchange-correlation 
potential is a rigid shift and the exact Vxc{t) is easily constructed. We also employed the 
adiabatic local spin density approximation (ALSDA) with self-interaction correction (SIC) 



24| . Since no reliable correlation potential is available for a one- dimensional electron system 
we only consider exchange. The Li norm of the deviation between the exact and TDDFT 
density is about 0.2 for the highly correlated system of u = 0.25. The ground state Kohn- 
Sham equation generates a set of excited virtual Kohn-Sham orbitals \n). Figure 1 shows 
a comparison of exact occupation probabilities and those obtained from the approximate 
S-matrix (eq. ()0.7j) ) by projecting '^^^^^'^{t) = $(xi, t)$(x2, t) onto Kohn-Sham Slater 
determinants. Shown are the occupation of the ground state (figure a, projection onto 
the Kohn-Sham determinant |0,0)), the first excited state (figure b, projection onto |0, 1)) 
and the second excited state (figure c, projection onto 0,2)). The second excited state 
\Ncm. = '^,nrei = 0) iuvolvcs a configuration mixture of at least two Kohn-Sham Slater de- 
terminants to be well represented (configurations |0,2) and |1,1)). Neither the projection 



8 



onto a single Kohn-Sham configuration state |0,2) (figure Ic) nor the projection onto the 
exact excited state (not shown) yields satisfactory transition probabilities. After the switch- 
off of the laser, the density undergoes oscillations resulting in time-dependent Hartree and 
exchange correlation potentials which give rise to oscillations in the occupation probabili- 
ties. These are the signatures of the "spurious cross-channel correlations" well-known from 
TDHF [l^ . For the present system, exact excited final states can be easily calculated. We 
have checked that cross-channel correlations persist when we project onto exact exit-channel 
states rather than Kohn-Sham determinants. Moreover, we also tested channel states ob- 
tained from a TDDFT linear-response equation. In this equation the exchange-correlation 
kernel of TDDFT was approximated by the exchange-only time-dependent optimized ef- 
fective potential 0, and the exact Kohn-Sham orbitals obtained by the exact ground 
state exchange-correlation potential have been used. The obtained excited-state wavefunc- 
tions, however, do not significantly differ from the single Kohn-Sham Slater determinants 
although the excitation energies are considerably improved compared to the Kohn-Sham en- 
ergy differences. With presently available approximations to the exchange- correlation kernel, 
TDDFT linear-response theory does not provide improved excited-state wavefunctions. We 
thus conclude that the projection amplitude of eq. ()0.7|) is not well-suited to determine an 
approximate S-matrix within TDDFT. 

To test the newly proposed read-out functional which depends only on the density, we 
use the exact time-dependent density n{x,t). In this way, errors due to the approximate 
exchange-correlation potential can be ruled out and the quality of the proposed functional 
for state-to-state transition probabilities can be directly assessed. The sum in eqs. ()().12j] . 
()0.13|) and ()0.14|) is truncated after the second excited state. Figure |21 shows a comparison 
of transition probabilities obtained by projecting the wavefunction according to eq. ()0.3|) 
(solid line) and by the new density functional of eq. ()0.14|) (dashed line) for the three lowest 
lying states. In the limit of t — > oo, i.e. as the averaging interval in eq. ()0.12|1 increases, 
the transition probability converges within the numerical accuracy towards the exact result. 
Numerical errors are due to the truncation of the sum over final states in eq. ()().14j) . Note 
that simply averaging over the cross-channel correlation in eq. ()0.7|) fsee figure^ would lead, 
in general, to incorrect results. Similar good agreement can be found for other initial states 
and other systems 26^ . 

In conclusion, we have investigated two different methods to extract state-to-state transi- 



9 



tion amplitudes from TDDFT calculations. In the first method, the correlated many-body 
wavef unction is approximated by a Slater determinant of the time-dependent Kohn-Sham 
orbitals. This approximate wave-function is projected onto appropriate final states (exact 
states or Kohn-Sham configuration states). The resulting state-to-state transition prob- 
abilities suffer oscillations after the switch-off of the external perturbation. The second 
read-out functional to calculate state-to-state transition probabilities directly involves the 
time-dependent densities and represents thus a well-suited density functional within the 
framework of TDDFT. The problem of cross-channel correlations can be avoided and well- 
defined transition probabilities can be determined in the asymptotic limit t — > oo. First 
results obtained by evaluating the read-out functional with densities resulting from the 
exact solution of a one-dimensional two-electron Schrodinger equation are very promising. 
The application of the read-out functional to double excitation of helium are currently being 
investigated. Work supported by FWF-SFB 016. 
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